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Abstract 

Exact solutions of the Einstein's field equations describing a spheri- 
cally symmetric cosmological model without a big bang or any other kind 
of singularity recently obtained by Dadhich and Patel (2000) are revis- 
ited. The matter content of the model is a shear-free perfect fluid with 
isotropic pressure and a radial heat flux. Three different exact solutions 
are obtained for both perfect fluid and fluid with bulk viscosity. It turns 
out that the cosmological rerm A{t) is a decreasing function of time, which 
is consistent with recent observations of type la supernovae. 

PACS number: 98.80.Es, 98.80.-k 

Keywords : non-singular models, cosmology, variable cosmological constant, 
causality principle 



^Corresponding Author 



1 



1 Introduction 



The problem of cosmological singularity is one of the most fundamental issues in 
modern theoretical cosmology. Due to the powerful singularity theorem [ll [2] , 
it was widely believed that cosmological models must have initial singularity. 
However, in 1990 Senovilla :3 obtained the first singularity-free cosmological 
perfect-fluid (with a realistic equation of state 3p — p) solution of the Einstein 
equation and since then the possibility of constructing regular cosmologies was 
renewed. The interest for regular cosmologies had stifled for nearly 30 years 
due to the powerful singularity theorems, which seemed to preclude such space- 
times under very general requirements, such as chronology protecting, energy 
and generic conditions. The open way to regular cosmologies was found in the 
violation of some technical premises of the theorems. The remarkable feature is 
the absence of an initial singularity, the curvature and matter invariants being 
regular and smooth everywhere. This corresponds to a cylindrically-symmetric 
spacetime filled with an isotropic radiation perfect fluid. For instance, it was 
shown by Chinea et al. [4 that the Senovilla spacetime did not possess a compact 
achronal set without edge and could not have closed trapped surfaces. However, 
the flrst results were not encouraging. The extension of the Senovilla solution 
to a family of spacetime left the set of regular models limited to a zero-measure 
subset surrounded by spacetime with Ricci and Weyl curvature singularities [5] . 
A thorough discussion of the model such type can be found in Senovilla ^ . This 
family is shown to be included in a wider class of separable cosmological models, 
which comprises FLRW universe [J. Other properties of these solutions, such 
as their inflationary behaviour, generalized Hubble law and the feasibility of 
constructing a realistic non-singular cosmological model are studied therein. 

A large family of non-singular cosmological models and generalization thereof 
have been considered but they all are cylindrically symmetric ^ — [10] . For 
practical cosmology the spherical symmetry, however, is more appropriate. It 
is therefore pertinent to seek spherically symmetric nonsingular models. The 
first model of this kind was obtained by Dadhich [11] with an imperfect fluid 
with a heat flux. The model satisfied all energy conditions and had no singu- 
larity of any kind. Dadhich et aL[12j also obtained a non-singular model with 
null radiation flux. These models are both inhomogeneous and anisotropic and 
have a typical behaviour beginning with two density at t — > —oo, contracting to 
high density at t — and then again expanding to low density at t —^ oo. An 
interesting feature of the spacetime metric of these models is that it contains an 
arbitrary function of time which can be constrained to comply with the demand 
of non-singularity and energy conditions, Dadhich and Raychaudhuri [13^ later 
showed how a particle choice of this function leads to a model of an ever exist- 
ing spherically symmetric universe, oscillating between two regular states, which 
involves blue shifts as in the quasi steady state cosmological model of Hoyle, 
Burbige and Narlikar |14j and is flUcd with a non-adiabatic fluid with anisotropic 
pressure and radial heat flux. These observations led to the search of spheri- 
cally symmetric singularity-free cosmological models with a perfect fluid source 
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characterized by isotropic pressure. Due to this search Tikekar [TS] constructed 
two spherically symmetric singularity-free relativistic cosmological models, de- 
scribing universes filled with non-adiabatic perfect fluid accompanied by heat 
flow along radial direction. Recently, many researchers [TB] — [lOj have studied 
non-singular cosmological models in different context. From a purely theoret- 
ical point of view, the investigation of nonsingular cosmological models gives 
invaluable insight into the spacetime structure, the inherent nonlinear character 
of gravity and its interaction with matter fields. As a by-product it also deep- 
ens our understanding of the singularity theorem, in particular the assumptions 
lying in their base [7]. 

Models with a dynamic cosmological term A(t) are becoming popular as they 
solve the cosmological-constant problem in a natural way. There are significant 
observational evidence for the detection of Einstein's cosmological constant, A 
or a component of material content of the universe, that varies slowly with time 
and space and so acts like A. New observations of type la supernovae (Garnavich 
et al. [H], Perlmutter et al. [22], Riess et al. (23], Schmidt et al. [IJ), cosmic 
microwave background (CMB) anisotropics (e.g., Lineweaver [25]), galaxy sur- 
veys (e.g., Lahav and Bridle [2S]), gravitational lensing (e.g., Chiba and Yoshii 
[27]), and the Lya forest (e.g., Weinberg et al. [2S]) argue for a nonzero cos- 
mological "constant" with riA(= K/iH"^) « 0.6 — 0.7. This quantity may not 
he constant as has been appreciated for 30 years (Bergmann [2S] , Wagoner [3D] , 
Linde [31], Kazanas [32]). Scalar fields (Dolgov [33], Abbott [34], Barr [35], Fee- 
ble and Ratra [36], Friemann et al. [37], Moffat [38|, Starobinsky [S^), tensor 
fields (Hawking [40|, Dolgov [41]), D-branes (Ellis, Mavromatos and Nanopoulos 
[42]). nonlocal effects (Banks [43], Linde [44]), wormholes (Coleman 05]), infla- 
tionary mechanisms (Brandenberger and Zhitnitsky [JS], Peebles and Vilenkin 
[47]). and cosmological perturbations (Abramo, Brandenberger and Mukhanov 
[48] ) have all been shown to give rise an effective cosmological term that de- 
cays with time. Earlier researches on this topic, are contained in Zeldovich [35], 
Weinberg [SD], Dolgov [5T], Bertolami [S2], Felten and Isaacman [S3], Charlton 
and Turner [5l], Sandaga [5S], Carroll, Press and Turner [SB]. Some of the re- 
cent discussions on the cosmological-constant "problem" and consequences on 
cosmology with a time-varying cosmological-constant have been discussed by 
Dolgov and Silk ^J\, Sahni and Starobinsky [58j, Peebles [5^, Padmanabhan 
[60] . Carroll [61], Vishwakarma [62], and Pradhan et aL[63]. This motivates us 
to study the cosmological models, where A varies with time. 

Recently Dadhich and Patel [64; obtained a shear-free nonsingular spherical 
model with heat flux. This model satisfies the weak and strong energy conditions 
and also has a physically acceptable fall-off behaviour in both r and t for physical 
and kinematic parameters. In this paper, motivated by the situation discussed 
above, we shall focus on the problem with varying cosmological constant in 
presence of perfect fluid and also in presence of bulk viscous fluid. We do this by 
extending the work of Dadhich and Patel [64] by including varying cosmological 
constant. The remainder of this paper is organized as follows. In Section 2 
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we give a description of the cosmological models with its dynamical equations 
and solve them under the initial conditions inspired by Dadhich and Patel. We 
also investigate three different cosmological models to different values of the 
function P{t) and discuss results for these regimes. Section 3 comprises bulk 
viscous universe. We present our discussions and conclusions in Section 4. 



2 A Perfect Fluid Universe Revisited 



In this section, we review the solutions obtained by Dadhich and Patel [64] . The 
metric of the model is given in the form 



ds^ = (r^ + Pf"dt^ - (r^ + P)2™ [dr^ + r^{de^ + sin^ Odcf)^)] , (1) 



where 



in particular, 



2to = 1 



< 0, 2n 



(2) 



Here P = P{t) which can be chosen freely. The Einstein field equations for a 
perfect fluid with time-dependent cosmological constant and a redial heat flux 
read: 



Rik - -j^Rgik + A^ifc = - 



[p + p)uiUk - pgik + 2 ili'^k + qkUi) 



(3) 



where we have set = 1, m^u* = 1 = —qiq^, QiU^ = 0, p and p denote the 
fluid density and isotropic pressure, and Qi is the radial heat flux vector. 
From Equations (H]) and ^ we obtain 



3m2p2 



(^j,2 _|_ p^2n+2 

m 

4 



4m{3P +{m+ l)?'^} 

fj.2 ^_ p\2m,+2 



A, 



^^2 _|_ p-j2m+2 



2{r^ + P)P + (3m -2n- 2)P^ 
[(m + n)P + n^r^] - A, 



4m(n + l)rP 



1 (^2 _^ p)n+2 ' 

where qi = qg}. The expansion and acceleration are obtained as 

{j.2 _|_ p)n+l ' 



(4) 

(5) 
(6) 

(7) 
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- (8) 

We have freedom of choosing the function P{t) so that to give a non-singular 
behaviour to the above parameters. As a matter of fact, there are muhiple 
choices (see, Dadhich and Patel |63]), for instant, P{t) = + b'^t'^, c? + e~^* , 
0^+52 cos Lot^ a} > \? . For all these choices it is observed that all the physical and 
kinematic parameters remain regular and finite for the entire range of variables. 
Note that the model admits an interesting oscillating behaviour in time, with 
oscillations between two finite regular states. Oscillating nonsingular models 
are quite novel and interesting in their own accord. 

For complete determinacy of the system, we assume an equation of state of the 
form 

V = 1P, 0<7<1, (9) 

where 7 is a constant. 



2.1 Model 1: 

We set P(t) = + , c? > }? . In this case the matter density p, the fluid 
pressure p, the heat flux parameter q and kinematic parameter of expansion 9 
are found to have the following expressions: 

12toW 4m [3(a2 + 62^2) + + 1)^2] 

-\- A, 



(r2 + a2 + 62i2)2«+2 (^2 + a2 _^ 52^2)2m+2 

..... ^. , J.2 + (3m - 2n - 1)62^2] 4 ^ ^)(-„2 _^ ^2^2^ ^ ^2^2j 

P= (^2 ^ ^2 ^2i2)2«+2 + (^2 +a2 .^52^2)2m+2 ^' 

_ 8m(n + l)r62t 

^ = (r2 +a2 +62i2)n+2' 

^ (r2+a2 + 62t2)n+i- 
Equations (fHI)) and pT|) . with the use of ([5]), reduce to 

x2 + (3m - 2n - 1 + 3m7)62t21 



(1 + 7)A 



(11) 
(12) 

(13) 



(r2 + a2 + 62i2)2«+2 

4 [(m + n + 3m7)(a2 + 62^2) + {^2 _^ j.^^^.^.^ _j_ ^^^|^2j 

+ (^2 .^^2.^ ^2^2)2^+2 • (l'^) 

From above equations it is evident that the matter density is always and ev- 
erywhere positive while positivity of pressure is ensured, if 3m ~ 2n < 1. The 
heat flux parameter q > for t < 0, q = for t — and g < for t > 0. 
Equation (jl3p implies that the model describes an expanding universe for t < 
with q > and a contracting universe for t > for g < 0, the switching from 
contracting phase to phase of expansion occurring at t = 0. 
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Figure 1: Variation of A with time for 2.1 Model 1. The values of parameters 
are: m = 1, n = 1 + VTz /2, 7 = 0.5, a = 2, b = 1 and r = 1. 



From Equation we observe that the cosmological constant A is a de- 
creasing function of time (see Figured]). We also observe that the A approaches 
a small and positive value at late times which is supported by recent type la 
supernova observations [21] — 



2.2 Model 2: 

We set P{t) — + e^^*" , > 6^. In this case the matter density p, the fluid 
pressure p, the heat flux parameter q and kinematic parameter of expansion 9 
are found to have the following expressions: 



12m262t2g-2fct^ 



(r2 +a2 + e-f'*') 



4m 



3(a2 + e-" ) + {m + lY 



(r2 + a2 + e-''*') 



^)^2^2m+2 



+ A, (15) 



Ambe 



P = 



(r2 + a2)(26i2 _ 1) _ g"''* + (3m - 2n)bt^e 



(^j,2 -|- a2 ^ g-fct2-)2n+2 

(m + n){a^ + e"''*^) + n2r2 



(r2 + a2 + e^&*^)2ni+2 

8nibr{n + l)te"^*^ 
= "(r2 +a2 -|-e-''t^)"+2' 

6m6ie"^*' 



-A, 



(r2 +a2 + e-''*') 



(16) 
(17) 
(18) 
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Figure 2: Variation of A with time for 2.2 Model 2. The parameters are: m 
1, n 1 - \/T7/2, 7 = 0.5, r = 1, a 2 and b = 1. 



By using Equation ^ and ehminating p{t) between (|15p and ([TB)l . we obtain 



(1+7)A = - 



Ambe 



{3(1 + 7)m - 2n}6t2e-''« + (^2 + a^){2bt^ - 1) 



{(1 + 3j)m + n}(a2 + e"''*') + {m(m + 1)7 + n'^}r^ 



(19) 



_|_ q2 _|_ g-6t2-j2m+2 

From above equations it is observed that the matter density is always and 
everywhere positive while positivity of pressure is ensured if 3to — 2n > 0. Thus 
it is also observed that the requirements of weak and strong energy conditions 
are fulfilled throughout the spacetime of this model. The dominant energy 
condition, which requires p > p, cannot, however, be satisfied: it is clearly 
violated for large r. Thus this model satisfies weak and strong but not the 
dominant energy condition. 

From Equation (fT9|) , it is observed that the A first decreases, reaches to a 
negative value then increases and becomes a constant small positive value (see 
Figure [2]) . This could play the role of dark energy. 



2.3 Model 3: 

We set P{t) ~ a^+b^ cos ujt, > b^. In this case the matter density p, the fiuid 
pressure p, the heat fiux parameter q and kinematic parameter of expansion 9 
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are found to have the following expressions: 

Sm^b^'w^ sin^ out Am[3{a'^ + 6^ cos cut) + {m + l)r'^] 



(r2 + 0,2 + p ^i)2n+2 (^2 ^ ^,2 _^ 52 (,Qg t^i)2m+2 

mb'^Lj^ [2{r^ + + b'^ cos Lot) coswt — (3to — 2rt — 2) sin^ wt 
4[(to + n)(o2 + 52 cos wt) + n2r2] 



A, (20) 



(r2 +a2 + 62 coswt)2'"+2 



A, (21) 



4m{n + l)rb'^uj sin ujt . , 

'^^ "(r2+a2 + ^,2coscji)"+2' ^ 



362ti; sinwi 



(23) 



(r2 +a2 + 52 coswi)"+i' 
Equations and (pij) . with the use of give 

rn62tj2 |^{3(i _f_ &2^^^ _ 2 (n + l)}sin2 - 2(r2 + a2 + ^2 coswi) coswi] 



(1+7)A = - 



(r2+a2 + 62 coswi)'"+' 

4 [{(1 + 37)m + n}(a2 + 52 cos wi) + {m(TO + 1)7 + n2}r2] 
^ (r2 +a2 + 62 cosa;t)2™+2 ' ^^"^^ 

From the above equations it is observed that the matter density is always 
positive whereas the pressure is non-negative if Sm — 2n < 0. Thus it can be 
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seen that the requirements of weak and strong energy conditions are fulfilled 
throughout the spacetime of this model but not the dominant energy condition. 
From equation (j23p the expansion parameter indicates that the universe of this 
model in the phase of contraction for 2a7r < ust < {2a + l)7r where a takes on 
integer values only. During the phase of contracting q < while during the 
expansion phase q > while q vanishing when switching from contraction to 
expansion occurs vice-versa. 

From equation (P^ . it is observed that the A is oscillating due to properties 
of sinusoidal functions (see Figure [3]). It is also worth- noting that the average 
value (with respect to one period) of A is positive. Here we will have negative 
equation of state at late times required to support the current acceleration of 
universe. 



3 Bulk Viscous Universe 

The equations of bulk viscosity can be obtained from the general relativistic 
field equation, we replace the effective pressure [15] 

(25) 

where p is the pressure due to perfect fluid present, ^ is the coefficient of bulk 
viscosity and 9 is the expansion scalar. Thus, given ^{t) we can solve for cos- 
mological parameters. In most of the investigations, involving bulk viscosity, 
it is assumed to be a simple power function of the energy density (Pavon [55] . 
Maartens Zimdahl [HZ]) 

m = Co/, (26) 

where and k are constants. If fc = 1, equation (|26p may correspond to a 
radiative fluid (Weinberg [SD]). However, more realistic models (Santos[5S|) are 
based on k lying in range < fc < ^. 
Introducing and (1^ into (O, we obtain 



P - 



4 



62 



2(r^ + P)P + (3to - 2n - 2)P 



p\2m+2 



[(to + 7i)P + nV^] - A. (27) 



3.1 Model 1: 

We set P{t) = + b^t^ , a? > . In this case we consider two following cases. 
3.1.1 Case I : solution for ^ = 

When fc = 0, Equation reduces to ^ = (constant) and hence Equation 
(|27)) with the help ^ and dTOj) reduces to the form 

_ Qmb'^iot 4mb^ [(2n + l)b^f - (r^ + g^)] 

( +VP~ (^^2 ^ ^2 ^ pi-2yn + l + (V2 +a2+62t2)2n+2 
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4 [(2m - n)(a2 + + {m(m + 1) - n^jr^] 

Eliminating p(t) between Equations (fTU]) and (pS)) . we get 

6m&2^oi 4m62 ^ + + {3(1 + 7)™ - 2n - 1)62^2] 



(1 + 7)A 



(r2 + a2 + 62i2)n+l (^2 + a2 + ^2^2)2^+2 

4 [{(1 + 37)m + n}(a2 + 52^2^ ^ |^2 _^ ^(-^ _,_ i)^}r2] 



2m+2 



(29) 



(r2 + a2 + 62^2) 
3.1.2 Case II : solution for ^ — £,qP 

When k — 1, Equation ([26|) reduces to ^ = ^oP and hence Equation (f27|) with 
the help of Equations ([9]) and (fTO|) takes the form 



1+7 



(r2 +a2 + 62^2)"+! 



4m62{(2n+ l)62i2 _ (^2 +a2)| 
P = - 



(r2 +a2 +62^2)2.1+2 



4 [(n - 2m){a? + fo^i^) + r2{n2 - m(m + 1)}] 

+ (r2+a2 + 62t2)2m+2 ' 

Eliminating p{t) between Equations ([501) and PU)) . we get 

[(1 + 7)(r2 + a2 + 62i2)n+l _ Q^^mbh]K = 

3mbH^ 3(a2 + ^2^2) + + l)r2 



(30) 



24^0^ & t X + ^2 _^ 52i2)2„+2 (^2 + ^2 _^ 52^2)2m+2 

4m&2 [r2 + a2 + {3(1 + 7)™ - 2n - l}62i2] 

(r2 +a2 + 62i2)«+i + 

4 [{(1 + 37)™ + n}(a2 + 62^2) _^ (j^2 _ ^j2 j^)^2j 

(r2+a2+62i2)2m-«+l " (^1) 

From Equations and (|3ip . we observe that the A is a decreasing function 
of time (see Figure 2]) , and it approaches a small positive value which is similar 
to the previously discussed Model 1 (Section [2+]) . 



3.2 Model 2: 

We set P{t) = a^ + e^''*^, > 5^. In this case we consider two following cases. 
3.2.1 Case I : solution for = 

When fc = 0, Equation ([26| reduces to ^ = (constant) and hence Equation 
([27|) with the help of Q and (fT5| reduces to the form 

(r2 + a2 + e-^*')'' 



(l+7)p = - ,^, , , + 
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Figure 4: Variation of A with time for 3.1 Model 1, case I (lower panel) and 
II (upper panel). Here m = 1, n = 1 + \/l7/2, 7 — 0.5, r = 1, a = 2, b = 
1 and ^0 = 1- 



Amhe 



(2m - n){a'^ + e^''*^) + {m(m + 1) - n^jr^ 



Eliminating p{t) between Equations (fTS]) and ([5^ . we get 



(32) 



(H-7)A = -- 



Ambe 



-bt-' 



{3(1 + 7)m - 2n}6t2e-b*' + (r^ + a^){2bt^ - 1) - e"''*' 



(■^,2 + a2 ^ g-6t2-)2ri+2 

{(1 + 37)m + n}{a'^ + er^^^) + {m{m + 1)7 + n^}r^ 



(r2 +a2 + e-''*') 
3.2.2 Case II : solution for ^ = ^oP 



bt^\2m+2 



(33) 



When fc = 1, Equation ((26)) reduces to f = ^oP and hence Equation ([271) takes 
the form 



7- 



(r^ + a2 + e"''* ) 
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Figure 5: Variation of A with time for 3.2 Model 2, cases I (lower panel) and 
II (upper panel). Here m = 1, n = 1 — \/l7/2, 7 — 0.5, r = 1, a = 2, b = 
1 and ^0 = 1- 



Ambe 



(H + a'^){2bt^ - 1) - e^"* + (3m - 2n)bt''e 



+ 



(r2 +a2 + e-''*')2»+2 
(m + n) [a? + e~^^^ ) + n^r 



Eliminating p{t) between Equations ([34]) and (fTS)) . we get 

[(1 + 7)(r2 + + e-''*')"+i - 6m6eoe"''*']A = 

3(a2 + e-''*') + (m + l)r2 



Ambe~''^^ 



Smb'^fe-^''*" 



(^2 _|_ q2 _|_ g-6t2'j2n+2 J^^2 _|_ q2 _|_ g-bt^ -12771+2 

{3(1 + 7)777 - 2n}bee-''*' + {r^ + a2)(26t2 _ i) _ e^''*' 



(^2 ^a2_^g-fct2)n+l 



{(1 + 37)777 + 77}(a2 + e"''* ) + {777(777 + 1)7 + 772}r2 



(r2 + a2 + g-6t2^2m-n+l 



(34) 



(35) 



From Equation (|33p we observe that the value of A first increases slowly and 
suddenly reaches to peak, then it has sharp decrease to a negative value, again, 
it has a slow increment and finally becomes a small positive constant value (see 
Figure [5] upper panel). From Equation (|35|) we observe that the A first decreases 
and then increases and finally approaches to a small positive constant. This 
could play the role of dark energy. 
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3.3 Model 3: 

We set P{t) = + 6^ cosujt, > 6^. In this case we consider two following 
cases. 

3.3.1 Case I : solution for ^ = 

When fc = 0, Equation ((26| reduces to ^ = (constant) and hence Equation 
((27)) with the help of Q and reduces to the form 

(1 + 7)P = 



(r2 +a2 + 52 coswt)"+i 
m&^w^ [{3to(62 - 1) + 2(n + l)}sin^ + 2{r^ + + 6^ cos wt) cos wt] 



(r2 +a2 + 62 coswt)2»+2 

4[(a2 + 62 cosLjt)(2m - n) - n2r2] 
(r2 +a2 + 62 cosu;i)2'»+2 " 

Eliminating p{t) between Equations (|20)) and ([36)) . we get 

3a;62go sintjt 
^ ^ (r2 +a2 + 62coscjt)«+i 

m62cj2 [{3(1 + 627)771 - 2n} sin^ cut - 2(r2 + a2 + 62 costjt) cosut] 
(r2 + a2 + 62 coswt)2»+2 

4 [{(1 + 37)m + n}(a2 + 52 coso^i) + {m(TO + 1)7 + n2}r2] 
(r2 + a2 + 62 cosa;t)2'"+2 



(36) 



(37) 



3.3.2 Case II : solution for ^ — S,op 

When k — 1, Equation (pS)) reduces to ^ = ^op and hence Equation P7|) takes 
the form 

3^0^^'-^ sinujt 



7- 



P 



(r2 + a2 + 52 coswt)"+i 

mlP'uj'^ [2(r2 + a2 + 62 coswt) coswt — (3m — 2n — 2) sin^ ujt\ 
(r2 +a2 + 62 coswt)2"+2 

4[(m + n)(a2 + 62 coso^i) + n2r2] 
^ (r2 + a2 + 62 cosu;i)2™+2 ' 

Eliminating p{t) between Equations ([571) (P0)l . we obtain 

[(1 + 7)(r2 +a2 + 52 coswt)"+i +3^o6^t^sinwt] A = 



(38) 



+3^0^ ^'T-'^ sinwi 



3m6^ti;2 sin^ 



(r2 + a2 + 62 cos wt)2"+2 



4{3(a2 + 62 cosujt) + (m + l)r2} 
(r2 + a2 + 52 cosa;i)2™+2 
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Figure 6: Variation of A with time for 3.3 Model 3, case I (lower panel) and II 
(upper panel). Here m = 1, n = 1 + \/T7 /2, 7 — 0.5, r = 1, a = 2, b = 1, ^0 = 
1,0; = 7r/6. 



vnlp-ij^ [2(r2 + 0^ + 62 cost^t) coscji - {3to(1 - 76^) - 2(n + 1)} sin^ 
^ (r2 + a2 + fe2 cos wt)"+i 

4 [{(1 + 37)771 + n}(a2 + 6^ cos wt) + (m + + l)r2] 

(r2 +a2 +62 cosu;t)2™-«+i ' ^^^^ 

From Equations ([57]) and ([5^ we observe that the A oscillates with time due 
to properties of sinusoidal functions, present in these equations. The nature of 
these models are same as already discussed in Model 3 (Section 12. 3p . 



4 Conclusions 

We have obtained a new class of spherically symmetric inhomogeneous cosmo- 
logical models with a perfect fluid and also a bulk viscous fluid as the source 
of matter with a radial heat flux without a big bang or any other singularity. 
These are the shear-free nonsingular models. These are inhomogeneous and 
hence accelerating but not shearing. It is the heat flux that combines with pres- 
sure gradient to avoid singularity. From the point of view of realistic cosmology, 
these models share with the standard FRW model, spherical symmetry and the 
absence of shear. 



The cosmological constant is a parameter describing the energy density of 
the vacuum (empty space) , and a potentially important contribution to the dy- 
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namical history of the universe. The physical interpretation of the cosmological 
constant as vacuum energy is supported by the existence of the "zero point" 
energy predicted by quantum mechanics. In quantum mechanics, particle and 
antiparticle pairs are consistently being created out of the vacuum. Even though 
these particles exist for only a short amount of time before annihilating each 
other they do give the vacuum a non-zero potential energy. In general relativity, 
all forms of energy should gravitate, including the energy of vacuum, hence the 
cosmological constant. A negative cosmological constant adds to the attractive 
gravity of matter, therefore universes with a negative cosmological constant are 
invariably doomed to re-collapse [52]. A positive cosmological constant resists 
the attractive gravity of matter due to its negative pressure. For most universes, 
the positive cosmological constant eventually dominates over the attraction of 
matter and drives the universe to expand exponentially |70| . 

The cosmological constants in all models given in Sections 2.1 and 3.1 are 
decreasing functions of time and they all approach a small and positive value 
at late times which are supported by the results from recent type la supernova 
observations recently obtained by the High-z Supernova Team and Supernova 
Cosmological Project (Garnavich et al. Perlmutter et al. [52], Riess et al. 

Schmidt et al. [H]). Thus, with our approach, we obtain a physically rele- 
vant decay law for the cosmological term unlike other investigators where adhoc 
laws were used to arrive at a mathematical expressions for the decaying vacuum 
energy. Our derived models provide a good agreement with the observational 
results. We have derived value for the cosmological constant A and attempted 
to formulate a physical interpretation for it. 

This paper adds a novel family of shear-free models to Senovilla's first model 
|3] and a large family of cylindrical nonsingular models |5l [7l [H [10] and a large 
family of spherical nonsingular models [TTJ [TU [M] avoiding cosmic singularity 
in the absence of shear. 
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